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The dark matter may be a composite particle that is accessible via a weakly coupled portal. If 
these hidden-sector states are produced at the Large Hadron Collider (LHC), they would undergo 
a QCD-likc shower. This would result in a spray of stable invisible dark matter along with unstable 
states that decay back to the Standard Model. Such “semi-visible” jets arise, for example, when 
their production and decay are driven by a leptophobic Z' resonance; the resulting signature is 
characterized by significant missing energy aligned along the direction of one of the jets. These 
events are vetoed by the current suite of searches employed by the LHC, resulting in low acceptance. 
This Letter will demonstrate that the transverse mass—computed using the final-state jets and the 
missing energy—provides a powerful discriminator between the signal and the QCD background. 
Assuming that the Z' couples to the Standard Model quarks with the same strength as the Z°, the 
proposed search can discover (exclude) Z' masses up to 2.5 TeV (3.5 TeV) with 100 fb _1 of 14 TeV 
data at the LHC. 


The existence of dark matter provides one of the 
strongest motivations for physics beyond the Standard 
Model, and its discovery is one of the core missions for 
the Large Hadron Collider (LHC) program. Under the 
assumption that the dark-matter particle is neutral and 
stable, it escapes the detector and manifests as missing 
transverse energy (dr)- The LHC collaborations have 
developed a comprehensive search strategy to look for 
signals with significant JZp, accompanied by jets and/or 
leptons (see, e.g. [1] for a review). These searches are 
typically cast in terms of a Simplified Model [2] for su¬ 
persymmetry or an effective theory of dark-matter inter¬ 
actions [3, 4]. Yet if one relaxes the assumption that the 
dark sector is weakly coupled, a new class of dark-matter 
signatures emerge that evade this entire suite of analyses. 
Namely, it is possible that the dark matter has been lurk¬ 
ing undercover within hadronic jets. The purpose of this 
Letter is to propose a straightforward discovery strategy 
for these “semi-visible” jets. 

Semi-visible jets may occur if the dark matter is the 
stable (or meta-stable) remnant of a more complicated 
dark sector. The dynamics of non-trivial dark sectors 
have been explored in many contexts, e.g. [5-19]. In 
these models, the dark sector contains a dark-matter can¬ 
didate^) and possibly new force carriers and/or matter 
fields. Note that we are agnostic about how much of 
the cosmological relic density is accounted for by this 
dark-matter candidate. Messenger states that couple the 
dark sector to the Standard Model (SM) can exist. If the 
messenger is accessible at colliders, dark-sector states can 
be produced, leading to unique signatures such as large 
particle multiplicities, displaced vertices, multiple reso¬ 
nances, and lepton or photon jets [20-33]. 

Another possibility is that the final state resulting from 
strongly coupled hidden sectors may contain a new type 
of jet object—a semi-visible jet. In this case, the dark 
matter is produced in a QCD-like parton shower along 
with other light degrees of freedom that decay hadroni- 
cally. The result is a multijet+^r signature where one 


of the jets is closely aligned with the JEIt- A corner¬ 
stone of the standard multijet+^T searches is to require 
a minimum angular separation between the jets and 
to remove QCD background contamination arising from 
jet-energy mis-measurement [34, 35]. This implies that 
events containing semi-visible jets have a low acceptance 
for the currently implemented suite of searches. 

To further illustrate this point, Fig. 1 compares se¬ 
lected observables for QCD with those for example 
weakly coupled and strongly coupled dark-matter mod¬ 
els. The weakly coupled model is derived from super- 
symmetric theories and results from pair production of 
1.5 TeV scalar quark partners. Each squark decays to a 
jet and 1 GeV neutral dark-matter particle. The signal 
from the strongly coupled model, which will be described 
more fully later, comes from the production of a 3 TeV 
resonance which then decays to a pair of dark-sector par¬ 
ticles that subsequently shower and hadronize, yielding 
semi-visible jets. Both these examples yield topologies 
with jets and missing energy. As the left panel shows, 
the weakly coupled (labeled WIMP) and strongly cou¬ 
pled (labeled semi-visible jet) dark-sector models produce 
considerable with tails that extend beyond the QCD 
distribution. However, A <f> = min {A^q^, A</q 2 ^ T }, 
where j-[ .2 are the two hardest jets, is different between 
these models, as illustrated in the right panel. The 
A <f> distribution falls relatively steeply for the strongly 
coupled case, while it remains relatively flat for the 
weakly coupled scenario. Typical LHC searches require 
A</> > 0.4 [34, 35]. After requiring ffp > 500 GeV and 
Acf > 0.4, the acceptance of the WIMP (semi-visible) 
example is ^70% (7%). We also verified that the ap un¬ 
derperforms on semi-visible jet signals; a standard cut of 
ap > 0.55 gives ~ 20% (3%) efficiency for the WIMP 
(semi-visible) scenario [36]. Razor analyses [37] may 
prove useful, but require optimization for senri-visible 
jets—see Supplementary Material. 

To regain sensitivity to final states containing semi- 
visible jets, the cut on the angular separation A<f> must be 
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FIG. 1: (left) The distribution of transverse missing energy Zt for the QCD background (solid blue), as well as the semi-visible 
jet (dashed red) and WIMP (dotted green) examples, (right) The distribution of A tj> = min {A(f>j 1 p T , Ar/>j 2 p T }, where j 1,2 are 
the two hardest jets. 


removed. This comes at the expense of an unsuppressed 
QCD multijet background, which must be eliminated us¬ 
ing other techniques. In this Letter, we focus on the case 
where the dark sector is accessed via a heavy resonance. 
In such scenarios, one can take advantage of structure 
in the transverse mass—calculated using the final-state 
jets and fZ ?—to distinguish the signal from QCD. The 
strategy employed is similar to others proposed for semi- 
visible Higgs decays [38]. 

We now introduce an example Hidden Valley [7] model 
that will enable us to analyze the LHC sensitivity for 
semi-visible jets. This model is presented for illustra¬ 
tion and concreteness; semi-visible jets will be among 
the LHC signatures for a vast class of dark-sector theo¬ 
ries. The messenger sector is described by a simple phe¬ 
nomenological model for a TeV-scale U( 1 )' gauge boson. 
The new leptophobic Z' gauge boson couples to the SM 
baryon current Jg M : 

£ 3 _ 1 Z'M* Z’^ _ I Mf, Zf Z'» - g%Y ^ JsV (1) 

Note that the Z' is treated as a Stueckelberg field— 
the Higgs sector has been neglected as it is not relevant 
for the LHC phenomenology discussed below; the addi¬ 
tional matter needed to render the t/(l) of baryon num¬ 
ber anomaly free is also ignored. 

The dark sector is an SU(2)d gauge theory with cou¬ 
pling ad and two fermionic quark flavors \i = X' 1,2 with 
masses M t . The dark quark coupling to the Z' is g^,. 
In general, the couplings g^, and g^ff 1 do not have to be 
comparable; we focus on the case where g^, is large so 
that the Z' decays frequently to the dark sector. 

The SU (2 ) d confines at a scale A d <C M z > ■ A QCD-like 
dark shower occurs when Mf ~ A 2 d so that many dark 
gluons and quarks are produced, which subsequently 
hadronize. Some of these dark hadrons are stable, while 
others decay back to the SM via an off-shell Z' . The 
detailed spectrum of the dark hadrons depends on non- 


perturbative physics. Nonetheless, some properties of 
the low-energy states can be inferred from symmetry 
arguments. There are two accidental symmetries: a 
dark-isospin number U( 1 ) 1-2 and a dark-baryon num¬ 
ber U( 1 ) 1 + 2 , where “1” and “2” refer to the x% flavor 
index. For example, the mesons X 1 X 1 and X 2 X 2 are not 
charged under either of these symmetries, and are thus 
unstable. The other mesons (xiX 2 > X 1 X 2 ) and baryons 
(XiX' 2 , X 1 X 2 ) are charged under U( l)i_ 2 and U( 1 ) 1 + 2 , 
respectively, and are stable. 

The spin of the dark mesons is also important. Sim¬ 
ilar to the g b (g c ) and T (J/ip) of the bottom (charm) 
system, the pseudoscalar and vector mesons should be 
degenerate. By naive degree-of-freedom counting, the 
vector mesons are produced roughly three times as often 
as the pseudoscalar mesons [29]. This impacts the phe¬ 
nomenology as the pseudoscalar decay is suppressed by 
a mass insertion, unlike the vector case. Therefore, the 
pseudoscalar decays are dominated by 6-quarks and are, 
in general, more displaced than the vector decays. The 
search strategy discussed below does not rely on 6-tags or 
displaced tracks; it may be possible to take advantage of 
the pseudoscalar meson decays with a more sophisticated 
analysis. 

Despite the myriad of possibilities for the dark sector, 
only certain parameters have a direct impact on the jet 
observables and missing transverse energy. The strength 
of the dark shower, parametrized by a d , plays a critical 
role. The coupling a d controls how many dark hadrons 
are emitted in the shower as well as their pt distribu¬ 
tions, which has a direct and measurable impact on the 
jet observables. In addition, the mass scale of the dark 
quarks affects the jet masses. 

The number of dark-matter particles produced in the 
shower impacts This effect can be parametrized as 

# of stable hadrons \ 

# of hadrons / 
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FIG. 2: Mass distributions after event selection cuts for the benchmark model in Table I, for various a d and ri lv . Mjj is the 
mass of the two large reclustered jets, Mt the transverse mass, and M mc the reconstructed Z' mass using all the dark-matter 
particles in the Monte Carlo. The r/ — rf> lego plots show the corresponding energy deposition in the detector. Red circles 
indicate visible SM hadrons, while the grey circles indicate undetected stable mesons. The crosses indicate the position of 
anti-^T R = 0.5 jets. The relative size of each circle and cross is set by the y / pr of the object. 


The value of r inv depends on the details of the dark-sector 
model. For the model described above with Mf = M| , 
the average proportion of the stable and unstable hadrons 
is equal, implying ri nv — 0.5. This assumes that the 
hadronization process is flavor-blind and that the dark 
quark masses are degenerate, and ignores baryon pro¬ 
duction, which is suppressed by a factor of 1/Afr, where 
N c is the number of dark colors. 

A mass splitting between the flavors can lead to vari¬ 
ations in r; nv . Assuming M 2 > Mi, in the Lund string 
model [39], fragmentation into heavier dark quark pairs 
is suppressed by the factor 


T = exp 


( 4tt| 

l ^ )■ 


( 3 ) 



DESCRIPTION 

BENCHMARK 

a x Br 

production rate 

80 fb 

M z , 

Z' pole mass 

3 TeV 

M d 

dark hadron mass scale 

20 GeV 

a d ( 1 TeV) 

running dark coupling 

0.2 

rinv 

fraction of stable hadrons 

0.3 


TABLE I: Parametrization for semi-visible jet search. 


Because of the exponential dependence of the fragmen¬ 
tation process, r; nv is very sensitive to small splittings of 
the mass parameters. As a result, fewer stable mesons 
are produced when M| — Mf > A. 2 d . This decreases the 
value of rinv below 0.5. To increase rj nv above 0.5, one 
can increase the number of flavors Nf, thereby enlarg¬ 
ing the number of stable mesons by Nf(Nf — 1), while 
only increasing the number of unstable mesons by Nf. 
Clearly, ri nv can take on any value between (0,1). 

Table I summarizes the five parameters that are most 
relevant for semi-visible jet observables. Three are sensi¬ 
tive to the details of the dark sector: the running dark- 
sector gauge coupling ad( 1 TeV), n nv , and the mass scale 
for the dark mesons M d . Note that by only including 
one value of Md, we are assuming that the LHC will 
be insensitive to the dark spectrum mass splittings, i.e., 
Mz> A d- Additionally, there are two portal parame¬ 
ters: the production rate er x Br and the Z' mass. 

To perform a detailed collider study, uu, dd —> Z' —» 
X ' X events were simulated for the 14 TeV LHC using 
PYTHIA8 [40] with the default CTEQ6 parton distribution 
functions. The dark-sector shower was simulated using 
the Hidden Valley Pythia module [28, 29], modified to 
include the running of ad as was done for [33]. Each 
meson had a probability ri nv to be a dark-matter par¬ 
ticle. The possible decays of dark baryons/mesons into 
each other were neglected. The resulting particles were 
processed through DELPHES3, with the CMS settings [41]. 

Anti-fcr R = 0.5 jets [42] were constructed and 
then reclustered into two large jets [43] using the Cam- 
bridge/Achen (CA) algorithm [44] with R = 1.1. In a 
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standard resonance search, one would use the invariant 
mass M?j = ( P 1 +P 2 ) 2 , where pi^ are the momenta of the 
two final large jets j 12 - However, the M n variable is not 
useful when there are a significant number of dark-matter 
particles in the shower. A variable that incorporates the 
missing momentum is the transverse mass: 

Mr = Ml + 2 ( v .\/;- y • ,4„ y.T - PTjj -Mr). (4) 

In a detector with perfect resolution, Mjj < Mt < Mz> ■ 
Figure 2 shows the distribution of Mjj , Mt and M mc af¬ 
ter event selection. M mc is the reconstructed Mz> com¬ 
puted from all the reclustered jets and truth-level dark- 
matter four-vectors. Mt in general yields a narrower, 
more prominent peak closer to M mc . The top panels of 
Fig. 2 show sample events for the different signals. The 
dark-sector particle multiplicity decreases for smaller ad- 
As r; nv is increased, the signal degrades because more 
stable mesons are produced and more information is lost. 

To estimate the reach at the LHC, we simulated 
60 x 10 6 QCD events, 5 x 10 6 W ± /Z + jj events, and 
5 x 10 6 tt events. All samples were binned in Ht in 
order to increase statistics in the high -Mt tails [45] us¬ 
ing Madgraph5 [46] at parton level and PYTHIA8 for the 
shower and hadronization. The dominant background af¬ 
ter event selection is QCD and W ± /Z+jj. For the signal, 
25000 events were generated for each choice of Mz' in in¬ 
crements of 500 GeV, using the benchmark parameters in 
Table I. An 8 TeV sample was used to validate the QCD 
background and limit-setting procedure [47] against the 
CMS dijet resonance search [48]. The JZt distribution 
was also validated [49]. 

Each event was required to have at least two R = 0.5 
anti -/ct jets with pt > 200 GeV and |r;| < 2.5, as well 
as JZt > 100 GeV. These pre-selection cuts model the 
impact of the trigger. Then, the following cut-flow was 
applied: 

• Recluster jets into R = 1.1 CA jets (j i, j 2 ); 

• Require \r]j 1 — r]j 2 \ < 1.1; 

• Require A</> < 1, where A (f> is the minimum az¬ 
imuthal angle between J^t and PTj 12 ; 

• Veto isolated e ±/with pt > 20 GeV, \r/\ < 2.4; 

• Require y t/Mt > 0.15. 

The R = 1.1 jets capture the wider radiation pattern 
expected from dark-shower dynamics. The cut on the 
pseudo-rapidity difference removes f-channel QCD [48, 
50] . The lepton veto and A (j> requirements suppress elec- 
troweak backgrounds. Finally, the y t/Mt cut effectively 
acts as a missing energy requirement; cutting on the di¬ 
mensionless ratio avoids sculpting the Mt distribution. 

After applying these cuts, a bump hunt was per¬ 
formed using Mt- Following the dijet resonance searches 
[48, 50, 51], the resulting background distribution was 


yfs = 14 TeV f Ldt = 100 fb 1 



FIG. 3: Expected a x Br for the signal benchmark in Table I. 
The shaded green (yellow) band corresponds to ±1(2) stan¬ 
dard deviations. The dashed black line shows the a x Br for 
a Z' with the same coupling to quarks as the SM Z°\ the 
dashed red line shows estimates 5<r discovery. The shaded 
purple region indicates where the vertices are displaced, as¬ 
suming g^i — 1 and that the shower is dominated by vector 
mesons. 


parametrized using a fitting function—see Supplemen¬ 
tary Material. Assuming the background exactly follows 
the fit obtained from simulation, the exclusion reach for 
the signal benchmark can be computed. Figure 3 shows 
the results for 100 fb -1 of 14 TeV LHC data as a function 
of M z > for the benchmark parameters (Table I). We as¬ 
sume a 10% width for the Z', as computed using the 
benchmark parameters. The production cross section 
times branching ratio for a Z 1 with the same coupling 
as the SM Z° is shown as a reference. A Z' with SM 
couplings can be discovered (excluded) up to masses of 
~ 2.5 TeV (3.5 TeV). 

We estimate that the dijet limit on cr x Br(Z' qq ) is 
comparable to the limit obtained for the dark-sector de¬ 
cay mode. For g z , ~ 1, the branching ratio to the dark 
sector varies from 80% to 50% along the expected exclu¬ 
sion bound as the Z' mass increases. Thus, the model 
would be discovered in the semi-visible jet channel before 
it would be observed by the irreducible dijet channel; this 
conclusion only gets stronger for more integrated lumi¬ 
nosity. 

We simulate prompt decays for the dark mesons. For a 
sufficiently heavy Z' and small couplings, the dark vector 
meson decays could yield displaced vertices. Requiring 
that the lab-frame decay length be < 0(1 mm), a lower 
bound on the couplings can be obtained: 


.4“ > io- 



fB f M z , 
V 10 \3 TeV 


/ 20 GeV\ 

J 


(5) 


where B ~ 10 is the average boost factor computed from 






5 


the benchmark simulation. Eq. (5) gives the lower purple 
region in Fig. 3. Modifications of the search strategy can 
still be effective in this region. 

This Letter proposed a new search strategy for the dis¬ 
covery of hidden-sector physics in resonance searches. In 
particular, the focus was on dark-sector showers that re¬ 
sult in novel semi-visible jets—objects that are composed 
of SM hadrons and dark matter. We argued that this 
generic signature could arise from a large class of strongly 
coupled dark-matter models. Furthermore, we gave a 
simplified parameterization that allows for a systematic 
treatment of the signature space. Finally, we provided 
expected exclusion limits using a bump hunt in trans¬ 
verse mass. A Z' with SM-size couplings to quarks could 
be discovered (excluded) up to ~ 2.5 TeV (3.5 TeV). 

There are two main extensions that can be explored. 
First, one can allow for leptons, photons, and/or heavy- 
flavor particles to be produced in the shower. Second, 
one can consider other production modes. In this case, 
the semi-visible jets may not be aligned with the 
and additional variables using jet substructure, along the 
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Supplemental Material 


A. Razor Analysis 


Inclusive shape analyses can provide an alternative approach in the hunt for semi-visible jets. One such example is 
the use of “razor” variables, which were originally designed to extract signals such as pair production of squarks that 
decay into jets and neutralinos [53]. The starting point is to recluster an event into two “mega-jets,” where the sum 
of the squared-invariant-mass of the jets is minimized [54]. Given these two mega-jets j\ , j- 2 , one can define the Mr 
variable: 

M R = \J(\Ph\ + \Ph\) 2 - (Pz,h +Pz,h) 2 ■ (SI) 


For general new-physics signals, Mr approximates the hard scale of the interaction. Additional discrimination between 
signal and background can be gained by considering the missing momentum JSr- For razor searches, the R 2 variable 


R 2 


1 

2 M| 


+ PTj 2 ) - Mt ■ (PT,j 1 


+ 


PT,h) 


(S2) 


is utilized. R 2 characterizes the separation between fZ? and jet momenta. 

Figure SI shows the distribution of (Mr, R 2 ) for the strongly coupled (labeled semi-visible jet) and weakly coupled 
(labeled WIMP) models in the main Letter (see Fig. 1). The WIMP signal is broadly distributed in the bulk of 
the (Mr,R 2 ) plane. On the other hand, semi-visible jets yield a distinctive edge at ( Mr,R 2 ) = (Mz>, 0), with a 
tail extending to R 2 > 0 and Mr < Mz> ■ Current razor analyses delineate exclusive signal regions in the (Mr,R 2 ) 
plane. For example, in the CMS dark-matter search that uses razor variables [37], the signal region is defined to 
have R 2 > 0.5 and Mr G (200, 300], (300,400], (400, 600] and (600, oo) GeV. Other analyses require less stringent cuts 
on R 2 ; for example, the supersymmetry analysis requires Mr > 400 GeV and 0.18 < R 2 < 0.5 for fully hadronic 
events [55]. None of these signal regions is ideal for the semi-visible jet signal, as they do not take advantage of the 
sharp kinematic feature in the Mr distribution. Re-optimizing the razor analyses may prove useful in distinguishing 
semi-visible jets from background; we leave this study for future work. 


Semi-visible Jets WIMP 



M r (TeV) M r (TeV) 


FIG. SI: The R 2 vs. Mr distribution for the strongly coupled (left) and weakly coupled (right) dark-matter models from 
Figure S2 in the main Letter. The color scale denotes the fraction of events in each bin. 
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B. Background Distributions 


Figure S2 shows the Mt distributions of the backgrounds before and after the final JZt/Mt cut. Requiring 
> 0.15 removes a significant fraction of the QCD background; afterwards, the W ± /Z°+jj and tt backgrounds 
dominate at large Mt- The background distribution is modeled using the fitting function 


f(x) = Po 


(1 _ x) v 1+P2 ln:c 

XP3+P4 In® 



(S3) 


where the pi are fit parameters; the best-fit curve is shown in dashed purple. The background distribution is smooth 
and the fit function provides a reasonable description in the region of interest. The Mt distribution for the benchmark 
signal listed in Table I is shown in solid black. The signal exhibits a prominent peak over the background and the 
f^T/M t cut significantly enhances the signal to background ratio. 



M t (GeV) M t (GeV) 

FIG. S2: Mt distribution before (left) and after (right) an J£t/Mt > 0.15 requirement. (All other selection criteria are applied 
for both panels.) The QCD, W ± /Z° + jj, and tt backgrounds are shown in red, green, and blue, respectively, and are stacked. 
The dashed purple curve is an analytic fit to the total background using Eq. S3. The solid black curve corresponds to the 
benchmark signal, with parameters listed in Table I. 








